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ABSTRACT 


The proposed Mars Science Laboratory (MSL) mission is intended to deliver a large rover to the Martian 
surface within 10 km of the target site. This paper presents computational fluid dynamics (CFD) predictions of 
forebody heating rates for two MSL entry configurations with fixed aerodynamic trim tabs. Results are compared to 
heating on a 10-deg sphere-cone reference geometry. All three heatshield geometries are designed to trim 
hypersonically at a 16 deg angle of attack in order to generate the lift-to-drag ratio (L/D) required for precision 
landing. Comparisons between CFD and tunnel data are generally in good agreement for each configuration, but the 
computations predict more flow separation and higher heating on a trim tab inclined 10 deg relative to the surface. 
CFD solutions at flight conditions were obtained using an 8-species Mars gas in chemical and thermal non- 
equilibrium. Laminar and Baldwin-Lomax solutions were used to estimate the effects of the trim tabs and 
turbulence on heating. A tab extending smoothly from the heatshield flank is not predicted to increase laminar or 
turbulent heating rates above the reference levels. Laminar heating on a tab deflected 10 deg from the conical 
heatshield is influenced by flow separation and is up to 35% above the baseline heating rate. The turbulent solution 
on the inclined tab configuration predicts attached flow and a 43% heating increase above the reference level. 


NOMENCLATURE 

A axial force (N) 

C D drag coefficient, D/ 'jy^VjSref 

CG center of gravity 

D heatshield diameter (m) or drag (N) 

D tah heatshield diameter with a tab ( m ) 
h altitude (km) 

L lift (N) 

L/D lift-to-drag ratio 

M Mach number 

MSL Mars Science Laboratory 

m aeroshell mass (kg) 

N normal force (N) 

p pressure (N/m 1 2 ) 

q heating rate ( W/cm 2 ) 

R„ nose radius (m) 

R s shoulder radius (m) 

Re Reynolds number 

S heatshield surface area (m 2 ) 

reference area, kD 2 / 4 (m 2 ) 

T temperature (K) 

TPS Thermal Protection System 
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t time from atmospheric interface ( s ) 

u velocity parallel to surface (m/s) 

V velocity (km/s) 

x/R b non-dimensional radial coordinate 

y + non-dimensional boundary layer coordinate 

a trim angle of attack (deg) 

5 boundary layer thickness (in) 

£ surface emissivity 

y flight path angle (deg) 

r\ normal distance to surface (m) 

At) grid cell height at wall (m) 

p viscosity (kg/m 2 -s) 

p density (kg/m 2 ) 

c 7 Stefan’s constant (5.67 x 10' 8 W/m 2 /K*) 

T shear stress (N/m 2 ) 

0 momentum thickness (m) 

8 C cone half angle (deg) 

Subscripts 

D heatshield diameter 

e boundary layer edge 

FR Fay-Riddell formula 

lam laminar 

turb turbulent 

w wall 

8 momentum thickness 

<» freestream 
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INTRODUCTION 


The Viking and Pathfinder missions successfully 
used ballistic entry trajectories to land on Mars with a 
target accuracy on the order of 100 km. Second 
generation missions, such as the proposed Mars Science 
Laboratory (MSL), will be designed to land with small 
footprints (< 10 km) at areas that are of particular 
scientific interest 1 ' 2 . The capability for precision 
landing, terrain hazard detection and avoidance, and 
soft touchdown proposed for MSL will allow targeting 
and exploration of sites that are otherwise unreachable. 

Guided atmospheric entry is proposed for MSL 
to achieve 10 km landing accuracy. Figure 1 depicts a 
sample entry, descent, and landing sequence using a 
two-parachute decelerator system. The guidance 
algorithm controls trajectory cross range and down 
range using vehicle lift in order to reach the desired 
supersonic parachute deployment conditions. Entry 
vehicle lift mitigates uncertainties in predicted entry 
states, atmospheric properties, aerodynamics, etc. that 
otherwise would result in a much larger landing 
footprint. A hypersonic lift-to-drag ratio (L/D) of 0.22- 
0.25 (trim angle of attack a - 16 deg) and ballistic 
coefficient (m/C D S ret ) less than 120 kg/m ' satisfy the 
delivery requirements for a 2005 launch. 


Atmospheric Entry (h = 125 km) 


**** Deploy Supersonic Parachute (M = 1 .6 - 2.28) 


i 



Deploy Subsonic Parachute (M = 0.8) 
& Jettison Heatshield 



Jettison Backshell & 
Supersonic Parachute 


Determine Altitude, 
Velocity, Hazards 


Powered Descent 

Hazard Detection & 
Avoidance 

^ Touchdown 


Figure 1. Sample MSL Entry, Descent, and Landing 
Sequence 


The baseline MSL aeroshell has a 4.05 -m 
diameter, 70 -deg sphere-cone heatshield that draws 
from Viking and Pathfinder heritage. Figure 2 shows a 
sample spacecraft cruise configuration with the 
aeroshell and cruise stage attached through the 
heatshield. The aeroshell will separate from the cruise 
stage prior to atmospheric interface and decelerate the 
spacecraft to supersonic speeds before parachute 
deployment. During the hypersonic portion of entry, 
the aeroshell must protect the spacecraft components 
from significant aerothermal loads. 



Figure 2. Sample MSL Cruise Configuration 


Blunt body lift is generated from axial force and 
is in the direction opposite the trim angle (Fig. 3). Trim 
conditions can be achieved with either radial center of 
gravity (CG) offset or an aerodynamic trim tab, or a 
combination of both. Using a fixed control flap may 
reduce the amount of ballast mass and relax internal 
packaging requirements needed to shift the radial CG. 



Figure 3. Blunt Body Trim Using a Fixed Tab 

If a trim tab is used instead of CG offset alone to 
achieve the required L/D , the flap heating environments 
are of particular interest. A trim tab would be exposed 
to high heating during entry, possesses a low thermal 
mass that can result in significant integrated heating, 
has a complex curved surface that presents 
manufacturing challenges, and may require adjoining 
dissimilar TPS materials that can recede unevenly. 
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Heatshield Configurations 

Aeroheating predictions are presented for three 
heatshield geometries: the reference axisymmetric 
shape and two configurations with fixed trim tabs (Fig. 
4). The flap on the Shelf configuration is a smooth 
extension of the forebody, whereas the Canted tab is 
inclined 10 deg relative to the flank. The Shelf and 
Canted geometries have 4.3% and 2.7% more surface 
area, respectively, than the Baseline shape. Both tabs 
were sized to eliminate the need for a CG offset. Table 
1 summarizes the geometric dimensions. 



Figure 4. Heatshield Configurations: Front (Top) 
and Side (Bottom) Views 


Table 1. Heatshield Dimensions 


Parameter 

Value 

e c 

70 

D 

4.05 

Dtab 

4.31 

Rn 

0.9854 

Rs 

0.0988 


This paper compares predicted aeroheating for 
the candidate heatshield geometries. Reference 3 
contains a detailed discussion of the Baseline 
aeroheating environments at multiple trajectory 
conditions. Select results are repeated here and 


compared to the Shelf and Canted heating. References 
4-9 present various analyses of the candidate MSL 
shapes: experimental aeroheating 3 4 ' 6 , supersonic 
aerodynamics 78 , and entry trajectory simulations 9 . 
CFD results in the present study are compared to data 
from Reference 6 for all three configurations. 

Atmospheric Entry Trajectory 

The MSL aeroheating analysis requires 
knowledge of the atmospheric entry path. The nominal, 
or expected, entry trajectory for the 2005 launch 
opportunity was calculated using the Program to 
Optimize Simulated Trajectories (POST) 10 . Figure 5 
shows the corresponding velocity-altitude map. 
Stagnation point heating rate, dynamic pressure, and 
Reynolds number are plotted versus time in Figure 6. 
The heating rate estimate is based on the Fay-Riddell 1 1 
formula for a sphere at a= 0. CFD solutions were 
obtained for the peak stagnation point heating time at 
the freestream conditions summarized in Table 2. 



Figure 5. Nominal Entry Trajectory for the 2005 
Launch Opportunity 



Figure 6. Stagnation Point Heating Rate, Dynamic 
Pressure, and Reynolds Number 
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Table 2. Peak Heating Point Freestream Conditions 


t (sec) 

103 

h (km) 

37.1 

V„ (m/s) 

4919 

P~ (kg/m 1 ) 

4.15 x 10' 4 * 

T m (K) 

159.3 

M m 

24.7 

Re D 

1.04 x 10 ft 


ANALYSIS 


Computational Approach 


Turbulence Modeling 

Turbulent LAURA solutions were obtained with 
the Baldwin-Lomax algebraic model 18,19 . The model is 
believed to give reasonable heating estimates for simple 
shapes, such as a sphere-cone. A transition location is 
required input and was specified to give conservative 
fully turbulent results. The Dhawan-Narashima 
bridging model 20 is employed in LAURA to blend the 
laminar and turbulent regions of the flowfield, but it is 
not intended to predict transition onset. More refined 
turbulence modeling (i. e. two-equation models) is 
recommended as the MSL aeroshell design matures, 
especially if a trim tab is used. 


CFD modeling will be combined with an 
experimental test program to predict the MSL flight 
vehicle aerothermal environments. Hypersonic flight 
conditions in a reacting Mars atmosphere gas cannot be 
recreated in an experimental environment, so the use of 
numerical tools is a critical part of the MSL aeroheating 
analysis. The continued improvement of grid 
generation methods, flow solvers, and computational 
resources makes it possible to solve three-dimensional, 
chemically-reacting flowfields that once were 
prohibitive. 

LAURA CFD Code 


Grid Adaptation 

LAURA is capable of grid alignment in which 
the computational mesh can be adapted to the boundary 
layer and bow shock based on user-defined parameters. 
Proper mesh resolution at the wall is necessary for 
reliable heating predictions. In LAURA, a user- 
specified cell Reynolds number controls the grid 
spacing at the wall: 


Re v .= 


paArj ^ 




(2) 


Aeroheating calculations were performed using 
the Langley Aerothermodynamic Upwind Relaxation 
Algorithm (LAURA) 12 CFD code. LAURA has been 
used previously to predict the aeroheating environments 
for other Mars applications 13,14 . LAURA models an 8- 
species Mars gas (C0 2 , CO, N 2 , 0 2 , NO, C, N, O) in 
chemical and thermal non-equilibrium using the Park- 
94 15 reaction rates. A finite-volume approach is used to 
solve the inviscid, thin-layer Navier-Stokes (TLNS), or 
full Navier-Stokes flowfield equations. The TLNS 
option was used for all calculations presented here. The 
code uses Roe’s averaging 16 for the inviscid fluxes with 
second-order corrections using Yee’s symmetric total 
variation diminishing (TVD) scheme 17 . 

A super-catalytic wall boundary condition was 
used in which the species mass fractions for C0 2 and 
N 2 are fixed at their freestream values of 0.97 and 0.03, 
respectively. A radiative-equilibrium temperature was 
imposed at the wall to satisfy the following relation: 

q w =£oT*, 

A fixed surface emissivity (e) of 0.9 was used for all 
solutions. The wall is assumed to radiate to a 
temperature of absolute zero. 


Experience has shown that reliable laminar heating 
predictions can be achieved with Re w = 10. Wall cell 
size for turbulent calculations was specified to insure 
sufficient resolution of the viscous sub-layer. A Re w 
value of 2 was used and generally resulted in v + of 
order 0. 1 : 


Singularity-free, structured volume grids were 
used for LAURA solutions on the Baseline, Shelf, and 
Canted configurations. Figure 7 shows surface and 
symmetry plane grids coarsened by a factor of two in 
each direction. Only half of the heatshield is modeled 
due to symmetry in the pitch plane. A 7-block mesh 
was used for the Baseline shape and 11-block grid 
systems were used for the Shelf and Canted 
configurations. The heatshield nose on each shape was 
mesned to avoid a singularity pole boundary that can 
introduce artificial discontinuities into the flowfield. 


+ 

y = 

n 




( 3 ) 


Computational Grids 
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Figure 7. Surface (Top) and Symmetry Plane 
(Bottom) Grids Used for CFD Solutions (Coarsened) 

The grids were built using GridGen 21 to 
construct the topology and surface distribution, 
3DGRAPE/AL 22 to generate the volume grid, and the 
Volume Grid Manipulator (VGM) 23 to enhance grid 
quality and impose boundary conditions along block 
interfaces. The Baseline surface grid has 3280 cells, 
and the Shelf and Canted surface grids each have 8048 
cells. The nominal meshes contain 64 cells in the 
surface normal direction. The grid distributions around 
the Shelf and Canted flaps were designed to accurately 
represent the surface geometry and sufficiently capture 
flowfield features caused by the tabs. 

RESULTS AND DISCUSSION 

LAURA perfect gas predictions are compared to 
select experimental data for the candidate MSL 
configurations. Laminar and turbulent predictions at 
flight conditions are shown for the peak heating time on 
the nominal trajectory. Select Baseline results are 


repeated from Reference 3 and compared with Shelf 
and Canted predictions to estimate trim tab effects. 

CFD vs. Tunnel Data 

Reference 6 reported experimental results from 
NASA LaRC’s 20-Inch Mach 6 Air Tunnel. Global 
heating data were obtained on the candidate geometries 
using the two-color, relative-intensity, phosphor 
thermography method. See Reference 6 for a detailed 
description of the facility, data acquisition, and test 
conditions. Figure 8 compares normalized heating data 
on the top half of the geometries at a = 20 deg and Re^ 
= 3.0 x 10 6 / ft. The Baseline and Shelf heating rates 
are essentially the same. Heating on the Canted model 
decreases ahead of the tab before increasing at the tab 
hinge line. It was inferred from this heating profile that 
the Canted tab caused flow separation. 



Figure 8. Comparison of Heating Rate Data at a = 
20 deg and Re „ = 3.0 x 10 6 / ft (from Ref. 6) 

LAURA solutions were run with perfect gas air 
and a constant wall temperature of 300 K. Comparisons 
between the computations and data from Figure 8 are 
shown in Figure 9. The estimated experimental 
uncertainty is +/-13% and is shown with the data. 
LAURA predictions are generally within the Baseline 
and Shelf data uncertainty. Predicted heating spikes at 
the shoulder do not show up in the data possibly 
because of model conduction effects. Computed 
heating on the Canted geometry agrees well with data 
on the conical flank, but tab heating is above the upper 
data boundary. The LAURA heating profile indicates a 
larger flow separation region ahead of the flap and 
stronger reattachment, which may explain the higher 
predicted tab heating levels. Further comparisons 
between CFD and tunnel data will be performed as the 
MSL aeroshell design matures. 
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Baseline 



Shelf 



Canted 



Figure 9. Comparisons of LAURA and Heating 
Rate Data at a = 20 deg and Re„ = 3.0 x 10 6 / ft 
(Data from Ref. 6) 


Laminar CFD Flight Predictions 

Computations at flight conditions were used to 
estimate the trim tab effects on the Shelf and Canted 
heating environments. Figure 10 compares LAURA 
laminar symmetry plane heating rates at the nominal 
trajectory peak heating time and a = 16 deg. The 
conical flank heating is generally higher on the 
windward side where the stagnation point is located. 
The nose experiences higher heating than the stagnation 
point even though they do not coincide. As expected, 
heating rates on the windward side are identical for all 
three geometries. Rapid flow expansion around the 
heatshield shoulder causes a spike in heating that can be 
mitigated by increasing the radius of curvature. 

The Shelf flap does not increase heating above 
the Baseline level. This is reasonable since the Shelf 
flap is an extension of the heatshield flank that does not 
introduce flow separation. The heating rate continues 
to decrease along the extra running length provided by 
the Shelf tab. Heating ahead of the Canted flap follows 
the Baseline level before passing through a pronounced 
drop and sharp rise along the flap corresponding to the 
starting and ending points of a separated flow region. 
The predicted heating rate on the Canted flap is as high 
as 49 W/cm\ which is approximately 35% higher than 
the laminar heating rate at the nose. 



Figure 10. Laminar Symmetry Plane Heating Rate 
Comparison (a = 16 deg) 

Figure 11 illustrates the beginning of flow 
separation ahead of the Canted flap and reattachment 
near the tab edge. The detached flow region spans the 
width of the flap and wraps around the tab sides as the 
flow expands. The recirculating flow pocket is 
centered at the tab hinge line and is followed by a weak 
shock at the reattachment point. 
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Figure 11. Laminar Pressure and Streamlines for 
the Canted Configuration ( a = 16 deg) 


A plot of surface pressure, streamlines, and 
heating rate on the Canted geometry is shown in Figure 
12. The separated flow ahead of the flap is easily 
identified in the heating rate distribution. Reduced 
heating in front of the tab corresponds with the 
separated How region. Peak heating occurs after 
reattachment as the flow accelerates around the tab 
edge. Heating rate on the comer is about 55 W/cm 2 , or 
12% higher than the symmetry plane value; a larger 
radius of curvature would reduce the heating. 



pipy] nmi! nn qjwWiBniii'n 


0 02 0 21 0.39 0.58 0.77 0 95 


2.7 13.1 23 5 34.0 44 4 54.9 


Figure 12. Laminar Pressure and Streamlines (Left) 
and Heating Rate (Right) on the Canted 
Configuration ( a = 16 deg) 


If the tab experiences elevated heating 
throughout the trajectory, the integrated heat load will 
be of concern because the TPS and structure in that area 
must perform at high temperatures. The heating 
problem is compounded by the fact that the exposed tab 
has low thermal mass and limited conduction paths 
through which heat is carried away from the flap area. 
Consequently, thermal soak may cause tab temperatures 
to continue rising after the applied heat pulse has ended. 

Figure 13 compares symmetry plane pressure 
between the configurations. The pressure is highest on 
the windward conical Hank where the stagnation point 
is located, and is somewhat lower on the opposite side. 
The detached How ahead of the Canted tab causes an 
increase in pressure, followed by another pressure jump 
near the flap hinge line. Pressure on the Canted tab is 
lower than the stagnation level. The higher pressure on 
the Canted Hap allows it to generate the desired trim 
angle of attack with less surface area than the Shelf 
geometry. 



Figure 13. Laminar Symmetry Plane Pressure 
Comparison (a = 16 deg) 

Canted Grid Refinement Study 

LAURA heating predictions were run until a 
grid-converged solution was obtained i. e. the heating is 
essentially unchanged with additional grid adaptations. 
A grid-independent solution is also desired, which 
means that the heating predictions do not change 
significantly for a more refined surface and/or volume 
grid. It was shown in Reference 3 that 64 normal cells 
are sufficient for grid-independent Baseline heating 
solutions. The same is also believed to be true for the 
Shelf heatshield because the trim flap is a smooth 
extension of the conical flank and does not cause flow 
separation. Thus, heating solutions on Baseline and 
Shelf are shown only for grids with 64 normal cells. 
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The Canted shape is unique because the trim tab 
induces flow separation whose recirculation pockets 
can be difficult to solve in a steady state solution. 
Consequently, the Canted was the most challenging of 
the three geometries for obtaining grid-converged and 
grid-independent laminar heating. In fact, global time- 
stepping was used at times to advance the flowfield 
solutions and allow the recirculating flow to stabilize. 

A grid refinement study was conducted to 
determine whether the reference Canted mesh with 64 
normal cells was sufficient to obtain reliable laminar 
heating predictions. The original surface grid was used 
and the number of surface normal cells was increased to 
96. Values of Re w = 10 and Re w ~ 1 were used to 
determine the effects of additional grid cells in the 
boundary layer and finer grid spacing at the wall. 

Figure 14 compares laminar heating predictions 
on Canted grids with 64 (Re w = 10) and 96 ( Re w = 10 
and Re w = 1) surface normal cells. Figure 15 shows an 
enlarged view of the flap region. The heating rate on 
the conical flank is unaffected by the number of grid 
cells or wall spacing. All three solutions exhibit flow 
detachment before the tab hinge line and a heating 
spike near the edge. Minor heating differences exist 
between the grids in the detached flow region and along 
the tab face. Table 3 summarizes the peak heating rates 
along the symmetry plane. The range of heating rates 
on the tab edge is 45.6-48.6 W/cm 2 , for a variation of 
less than +/-3.5% from the mean value. This range of 
heating rates is considered to be within the CFD model 
uncertainties and small enough to satisfy grid 
independent heating with respect to surface normal grid 
spacing. Ideally, solutions on a finer surface mesh also 
should be performed, but time and resource constraints 
prevented such an analysis in the present study. 



Figure 14. Laminar Grid Refinement Study on the 
Canted Configuration 



Figure 15. Laminar Grid Refinement Study on the 
Canted Configuration (Enlarged View) 


Table 3. Peak Laminar Symmetry Plane Heating 
Rates on Canted Configuration Grids 


No. of Surface 
Normal Cells 

Re K 

Peak q w (W/cm 2 ) 

64 

10 

48.6 

96 

10 

45.6 

96 

1 

47.6 


Turbulent CFD Flight Predictions 

Reference 3 showed that turbulent transition is 
likely to occur before the peak heating point on the 
nominal trajectory. Natural transition using a 
momentum-thickness Reynolds number criteria of Re 0 
> 200 and transition due to heatshield penetrations were 
both predicted to result in a turbulent boundary layer. 
Re e is commonly used as an indicator of smooth body 
transition and is based on the boundary layer edge 
properties and momentum thickness (6): 


Re d = P^- 


1 - — 


9 = j 

0 Pe^e 


PV 


( 4 ) 

( 5 ) 


The boundary layer edge is defined as the location 
where total enthalpy is 99.5% of the ffeestream value. 


If turbulence does occur before the time of peak 
heating, it will significantly augment the heating rates 
and integrated heat load. An estimate of the turbulence 
effects was accomplished using Baldwin-Lomax 
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computaiions on all three configurations. Reference 3 
reported more than a 70% turbulent increase in heating 
above the laminar nose level on the Baseline geometry. 
Those results are repeated here and compared to Shelf 
and Canted predictions. 

T urbulent heating rale along the symmetry plane 
of each configuration is shown in Figure 16. Contrary 
to the laminar profiles, the turbulent heating is higher 
on the leeward flank where boundary layer running 
length is largest. Again, the Baseline and Shelf heating 
rates are similar and peak at about 63 W/cm 2 . The 
longer running length on the Shelf heatshield does not 
result in additional turbulent heating augmentation. 
Turbulent heating on the Canted shape is identical to 
the Baseline levels on the commonly shared geometry. 
At the Canted tab hinge line, the heating immediately 
increases due to the surface inclination change. The 
heating rate across the entire tab face is elevated above 
the Baseline turbulent level. The heating rate on the 
Canted flap is 83-90 W/cm ", or 32-43% above the 
Baseline level and 69-84% higher than the Canted 
laminar peak value of 49 W/cm 2 . 



Figure 16. Turbulent Symmetry Plane Heating Rate 
Comparison (a = 16 deg) 

Elevated heating over the entire Canted flap is 
the result of attached turbulent flow (Fig. 17). The 
presence of a fuller and thicker turbulent boundary 
layer prevents the development of flow separation 
ahead of the Canted flap, contrary to the laminar case. 
This is a reasonable result since a turbulent boundary 
layer has more momentum (compared to a laminar one) 
with which to overcome a negative pressure gradient. 
The oblique shock emanating from the flap hinge line is 
much stronger than the laminar reattachment shock 
shown previously. 



Figure 17. Turbulent Pressure and Streamlines for 
the Canted Configuration (a = 16 deg) 


Surface pressure, streamlines, and heating rate 
on the Canted geometry in Figure 18 illustrate how 
attached turbulent flow corresponds to high heating 
rates on the tab face. The flap causes an immediate 
jump in pressure at the hinge line. The same is true of 
the heating rate distribution. Similar to the laminar 
case, the heating rate at the tab outer edge (98 W/cm 2 ) is 
slightly higher than the symmetry plane value and may 
be reduced with a larger radius of curvature. 



p/p V 2 IUI9DJLJH i| i wi i m 2 )M \\M 

0.02 0.25 0.47 0.70 0.92 1.15 5.6 24.1 42.6 61.0 79.5 97 . 


Figure 18. Turbulent Pressure and Streamlines 
(Left) and Heating Rate (Right) on the Canted 
Configuration (a = 16 deg) 


9 

American Institute of Aeronautics and Astronautics 



The turbulent surface pressures for each 
configuration are shown in Figure 19. Pressure 
distributions on the Baseline and Shelf geometries are 
similar to the laminar profiles. The Canted profile 
shows a much larger pressure jump with attached flow 
than was predicted for the detached laminar flow case. 
The pressure along the Canted tab is higher than the 
stagnation level. If the flow does remain attached on 
the Canted geometry, higher pressure on the tab would 
improve its aerodynamic effectiveness. 



Figure 19. Turbulent Symmetry Plane Pressure 
Comparison (a = 16 deg) 

Laminar and turbulent heating rates for the Shelf 
and Canted shapes are repeated in Figure 20. 
Turbulence is predicted to increase heating over the 
entire leeward flank for both heatshields. Windward 
side heating is predicted to increase slightly with 
turbulent conditions. Table 4 summarizes the peak 
symmetry plane laminar and turbulent heating rates for 
all configurations. Turbulent transition augments the 
Baseline and Shelf heating by the same magnitude of 
75% above the laminar nose level. The increase in 
peak heating on the Canted configuration is 84%. 

There is evidence that worst-case tab heating 
with flow separation may occur with a laminar 
incoming boundary layer followed by transitional 
reattachment 24 . It has been suggested that algebraic 
turbulence models can give reasonable heating 
estimates by specifying a transition location at 
reattachment. Furthermore, if the incoming flow is 
turbulent and is predicted to detach, more advanced 
turbulence models show promise in predicting heating 
rates on a deflected flap 25 . 




Figure 20. Comparisons of Laminar and Turbulent 
Heating Rates on the Shelf and Canted 
Configurations (a = 16 deg) 


Table 4. Peak Laminar and Turbulent Symmetry 
Plane Heating Rates (* denotes nose value) 


Configuration 

Q\am 

(W/cm 2 ) 

< 7 turb 

(W/cm 2 ) 

tflam turb 

Baseline 

36* 

63 

1.75 

Shelf 

36* 

63 

1.75 

Canted 

49 

90 

1.84 


The current analysis was used to estimate the 
effects of aerodynamic trim tabs on MSL heating. 
Numerical modeling of chemically reacting flowfields 
over deflected control flaps remains a challenging task, 
especially if flow separation exists. The prediction of 
attached turbulent flow for a modest control flap angle 
of 10 deg is plausible given that the turbulent boundary 
layer has more momentum to overcome a negative 
pressure gradient. If detached turbulent flow is 
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suspected to occur, more advanced techniques are 
recommended beyond the Baldwin-Lomax model. 
Additional testing and more refined models will be used 
to further analyze MSL turbulent heating, especially if a 
deflected flap is used to trim the entry vehicle. 

SUMMARY AND CONCLUSIONS 

Computational aeroheating predictions were 
obtained for candidate Mars Science Laboratory 
heatshields designed for a hypersonic lift-to-drag ratio 
of 0.22-0.25 at a trim angle of attack of 1 6 deg. The 
reference Baseline 70 -deg sphere-cone geometry 
achieves trim through an offset radial center of gravity. 
The Shelf and Canted geometries have fixed tabs that 
generate trim conditions through asymmetric loading. 
The Shelf tab extends smoothly from the heatshield 
flank and the Canted flap is deflected 10 deg from the 
conical heatshield surface. 

Perfect gas CFD solutions were compared to 
experimental results and show reasonably good 
agreement. CFD solutions on the Canted geometry 
predict more extensive flow separation and higher 
heating on the flap compared to data. Flight predictions 
for atmospheric entry from a 2005 launch were 
computed with the LAURA code using an 8-species 
Mars gas in thermo-chemical non-equilibrium. A 
radiative-equilibrium temperature boundary condition 
was specified at the wall. Laminar and Baldwin-Lomax 
turbulent solutions were obtained at the peak heating 
point on the nominal entry trajectory. 

Laminar heating on the Shelf configuration is 
identical to the Baseline level. Also, the additional 
running length on the Shelf geometry does not result in 
increased turbulent heating above the Baseline value of 
63 W/cm 2 . Peak heating rates on the Canted geometry 
are higher than the Baseline levels for both laminar and 
turbulent flow conditions. In the laminar case, the flap 
is predicted to cause flow separation ahead of the tab 
and reattachment on the flap surface. Peak heating 
rates on the flap edge are as high as 49 W/cm 2 , which is 
35% higher than the laminar nose heating rate. For 
turbulent conditions, the boundary layer is predicted to 
remain attached across the Canted tab hinge line. Peak 
heating rate is 43% higher than the Baseline level and 
84% above the Canted peak laminar heating rates. 

As the MSL aeroshell design matures, both 
experimental and computational methods will be used 
to predict aeroheating environments for TPS material 
selection and sizing. Comparisons between data and 
predictions will be used to validate the computational 
models, and CFD solutions at flight conditions will be 
used to predict worst-case heating levels. The 


predictions will include uncertainties in CFD modeling, 
entry states, atmospheric properties, tunnel-to-flight 
extrapolation, etc. for conservatism. The resulting TPS 
will be a robust design that has margin beyond the 
worst-case aeroheating environments. 
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